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ABSTRACT 

Extreme scattering events (ESEs) in the interstellar medium (ISM) were first observed in regular 
flux measurements of compact extragalactic sources. They are characterized by a flux variation over a 
period of weeks, suggesting the passage of a “diverging plasma lens” across the line of sight. Modeling 
the refraction of such a lens indicates that the structure size must be of order AU and the electron 
density of order 10s of cm“^. Similar structures have been observed in measurements of pulsar 
intensity scintillation and group delay. Here we report observations of two ESEs showing increases 
in both intensity scintillation and dispersion made with the Parkes Pulsar Timing Array (PPTA). 

These allow us to make more complete models of the ESE, including an estimate of the “outer-scale” 
of the turbulence in the plasma lens. These observations show clearly that the ESE structure is fully 
turbulent on an AU scale. They provide some support for the idea that the structures are extended 
along the line of sight, such as would be the case for a scattering shell. The dispersion measurements 
also show a variety of AU scale structures which would not be called ESEs, yet involve electron density 
variations typical of ESEs and likely have the same origin. 

Subject headings: pulsars: general — ISM: structure — methods: data analysis 


1. INTRODUCTION 

The original observation of ESEs (Fiedler et al, 1987), 
shown in Figure I, suggests the passage of a diverging 
lens across the line of sight, as the flux is weaker in the 
middle of the event and piles up towards the edge. A 
convex “blob” of high density plasma would provide such 
a lens. The refractive index variations would be smaller 
at a higher observing frequency and the observations are 
consistent with this behavior, so it has been assumed 
that this is the basic mechanism causing the ESE event. 

The observations attracted wide attention and many 
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reports of similar structures in the ISM have been pub¬ 
lished. These reports include pulsar observations of cor¬ 
related fluctuations in group delay and flux (Cognard 
et al, 1993, Lestrade et al, 1998), persistent phase gra¬ 
dients (Gupta et al, 1994), enhanced diffractive scintil¬ 
lation (Stinebring et al, 2001; Hill et al, 2005; Brisken 
et al., 2010); enhanced angular broadening (Lazio et al, 
2000) and increased dispersion (Keith et al, 2013). How¬ 
ever observations of both enhanced diffractive scattering 
and dispersion have not previously been reported for the 
same event. The combination allows one to make a more 
thorough analysis of the ESE including an estimate of 
the “outer scale” of the turbulence in the plasma. It has 
become common to use the term ESE for both the event 
and the plasma blob responsible for it. We will follow 
this convention. 

Pulsar timing arrays are designed to make precise mea¬ 
surements of the group delay of an array of pulsars every 
few weeks (see, for instance, Manchester et al. 2013 and 
references therein) for decades. The primary objective of 
PTAs is the direct detection of gravitational waves with 
periods of the order of a decade and for many pulsars 
the primary noise source is fluctuations in the electron 
density of the ISM. To correct for the ISM noise, PTAs 
also make precise measurements of the dispersion in the 
ISM every few weeks. These observations provide an ex¬ 
cellent window of opportunity to view fluctuations of the 
ISM on a scale of AU in general, and to study ESEs in 
particular. 

Pulsar observations are made with filter-banks because 
the change in dispersive delay over a typical observing 
band is much greater than the pulse width or even the 
pulse period, so the different frequency channels must be 
aligned before they can be averaged over the observing 
band. This provides a good measurement of the total 
electron column density, bnt for PTA purposes it must 
be supplemented by observations in two different bands 
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Figure 1. Original ESE observations of Fiedler et al (1987). 

separated more widely than any existing receiver band¬ 
width. Inter-band measurements of millisecond pulsars 
can achieve accuracies of 1:10® in dispersion. 

Filter-bank observations are also ideal for measuring 
the diffractive scintillation because the “dynamic spec¬ 
trum” of intensity scintillations is a two-dimensional ran¬ 
dom process with characteristic “scales” in both fre¬ 
quency, t'o, and time, tq. These scales depend on the 
strength of scattering, the distance to the pulsar, the 
velocity of the pulsar and the location of the scattering 
region. The distance and proper motion of the pulsar 
can often be obtained from the timing model. This is an 
astrometric model which is fitted to the time of arrival 
(ToA) of the pulses. For nearby pulsars the distance can 
be determined directly from the parallax, which causes 
a biannual sine wave in the ToAs. Otherwise the dis¬ 
tance can be inferred from the measured dispersion us¬ 
ing a model of the Galactic electron density (Taylor and 
Cordes, 1993, Cordes and Lazio, 2002). The location of 
the screen and the strength of scattering can then be 
estimated from the measured and tq. 

With measurements of the electron column density 
(DM) from the dispersion, and the level of scattering 
from z/Q, and tq, we can model the turbulence in the ob¬ 
served ESEs and compare this model with the “average” 
turbulence in the entire line of sight. In a general sense, 
it is surprising that one can detect structures in column 
density or scattering which are as small as an AU in a 
line of sight that is 100s of pc long. The exceptional ac¬ 
curacy of DM measurement provided by inter-band mea¬ 
surements make this possible. Scattering measurements 
cannot be made with accuracy greater than ^ 20%, but 
we will show that AU structures can cause exceptionally 
strong scattering if they are fully turbulent, i.e. if the 
fluctuations in electron density rie are comparable with 
the mean rie on an AU spatial scale. 

2. SCATTERING IN THE ISM 

A radio wave propagating through an ESE in the ISM 
will be scattered into an angular spectrum of plane waves 
B(0) by fluctuations in Ue transverse to the line of sight. 
The He fluctuations cause only phase changes in the ra¬ 
dio wave, but as the scattered waves propagate to the 
Earth they interfere, causing intensity scintillations to 
build up with distance from the ESE. When the scin¬ 
tillations are weak, i.e. the intensity fluctuations are 
small and the bandwidth is broad, the intensity fluc¬ 
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tuations will have a single characteristic scale which is 
of order of the Fresnel scale. As the scintillation gets 
stronger the characteristic scale bifurcates into a diffrac¬ 
tive scale, which decreases and a refractive scale which 
increases. The geometric mean of the two scales remains 
the same as the weak scattering scale. The diffractive 
scintillation causes a pulse delay Tdei which depends on 
6. This Tdei will de-correlate the intensity fluctuations 
at different wavelengths, giving a coherence bandwidth 
t'o = ll2'K{Tdei) ■ The refractive scintillations remain rel¬ 
atively broadly correlated over wavelength. Scintillations 
are normally observed as a time series, as the ISM drifts 
across the line of sight with a transverse velocity Ves- For 
the observations discussed here, and for typical pulsars 
at frequencies around 1 GHz, the diffractive scale is of 
order 10 min and the refractive scale is of order 1 day. In 
our observations, of the order of an hour, one sees diffrac¬ 
tive scintillations only, but the mean flux from one day 
to the next will vary due to the refractive scintillations. 
For a discussion of scattering see Coles et al. (2010) and 
for application to the ISM see Rickett (1990). 

2.1. Diffractive Intensity Scintillation 

Here we will refer to the scattering plasma as a gener¬ 
alized “screen”, which may be an ESE or simply a slab of 
ISM. When the pulsar is at a distance L and the screen 
is at distance (^L from the pulsar we can write 14ff as 

Kff = VeC + Vp{l — 0 — VisM- (1) 

Here all velocities are projected onto the celestial sphere. 
Tdei can be written as a function of 6 given the screen 
location. A wave scattered at angle 6 will leave the pulsar 
at angle 6p = {\ — (f)0 and will arrive at the Earth at 
angle 0e = C^- 

Tdei = C(1 - O0'^L/2c (2) 

If the phase fluctuations in the screen have stationary 
gaussian differences the autocovariance of the field at the 
screen Pe{i^ = (e(r')e*(r'-|-r)) is related to the structure 
function of the screen phase D^{r) = {{4>(r') — (j){r -|- 
r')Y) by /Oe(f0 = exp(—0.5D,^(r)). If the turbulence is 
Kolmogorov 

D 4 >{r) = (r/so)®/® for r < Sout and (3) 
= (so-ut/so)®'^^ for r > ^OUt ■ 

The brightness distribution B{0) is the Fourier trans¬ 
form of Peif), 

B{0) = jj pe{r) exp{—j2Trr- 0/\) cPr (4) 

Both Pe{i^ and B{0) are almost gaussian, so if we define 
the 1/^/e width of pe as sq then the I/\/e width of B{0) 
is 00 = A/27rso, and D^{so) = I. 

The bandwidth of the intensity fluctuations is then 
t'o = 1/2ttT dei i0o)- All the scintillation parameters, in¬ 
cluding pq and Tdei, are slowly varying functions of fre¬ 
quency. We can define them as approximately constant 
over a narrow band centered on pm ■ 

If vq < vm the scattering is strong and 

^'0 _ 2 / sq \ ^ 

VM 


(5) 
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where Vf = yj AL/27r is the Fresnel scale at vm- This is 
the case for all the observations in the PPTA except for 
those of PSR J0437-4715. 

The intensity time scale tq (at 1/e) is defined by sq = 
V/ff Tq , so with measurements of vq and tq and knowledge 
of L, one can find the location of the screen and the 
spatial scale sq. If we neglect Ve and Vjsm, which is 
often reasonable, we obtain 

C/(l - C) = {VpTolrff{2vMlj^o) and (6) 

So = Vp(l-C)ro. (7) 

If Ve is important one must add VeC/{1 — C) to Vp and 
the solution can be obtained iteratively. 

2.2. Outer Scale Model 

The outer scale of any turbulent system is difficult to 
measure but observations suggest that in most cases the 
outer scale is comparable with the dimensions of the sys¬ 
tem. Thus the outer scale Sout of an ESE should be com¬ 
parable with the smallest dimension of the ESE. The 
phase structure function limit for r > Sout is equal to 
twice the phase variance, D^{sout) ~ 2(^^) by defini¬ 
tion. If the path length Wz in the scattering structure 
equals Sout then because the rms density 

equals the mean density. This can be used to relate 
to {(j}) for thicker screens in which Wz = Ns out- In this 
case = {(f))'^/N. A spherical ESE would limit at 
Dcf,{Wz) = An ESE extended along the line of 

sight, such as a shell model, would have Sout ~ Wt the 
transverse width. The structure function would saturate 
at D^{Wt) = 2{4,)^{WtlWz). 

Observationally we can measure (</>) for the ESE 
through the change in group delay, i.e. (</>) = 2TTVM5tgrp 
and So from the intensity scintillation. We can also mea¬ 
sure Wt from the time it takes to cross the line of sight, 
but we don’t know Wz ■ liWt < Wz then Sout = Wt and 

{Wt/sof/^ = 2{2EVM5tgrpf{WtlWz). ( 8 ) 

We can solve this expression for Wz , which will be valid 
only if Wz > Wt ■ Otherwise one would have to set Sout = 
Wz and solve 

{Wz/st,f/^ = 2{2nvM5tgrpf ( 9 ) 

for Wz- The latter case did not occur in our observations. 

This analysis can also be applied to the entire line of 
sight, excluding the ESE. The scattering would be mod¬ 
eled as the superposition of N independent slabs of depth 
Sout distributed randomly over the line of sight. The C 
estimated would be a weighted average of the slab loca¬ 
tions. In this case we derive tgrp from the mean DM and 
solve 

{Sout/Sof/'^ = 2{2TTVMtgrpY{Sout/L). ( 10 ) 

for an average Sout- Such an analysis neglects a number 
of biassing effects and could be expected only to provide 
the correct order of magnitude for Sout- 

3. DATA ANALYSIS 

The dynamic spectra were measured with standard 
procedures in psrchive (Hotan, van Straten & Manch¬ 
ester 2004) which provide a flux estimate by fitting a tem¬ 
plate to the measured pulse profile. This is done for each 


frequency channel, m = I to M and each sub-integration 
n = I to N. The sub-integration time was I min and the 
observing time was usually 64 min. The auto-covariance 
C{t',p') of the dynamic spectrum is estimated 

using 

C{n ,m ) - EE S{n,m)S{n+n',m+m'). (II) 

This estimator multiplies the covariance by a triangle 
function but keeps the errors roughly independent of lag. 
We eliminate the bias caused by this triangle by includ¬ 
ing the triangle in the model fit as shown below. This 
temporal model is well supported theoretically, but the 
frequency model is only a rough approximation. A gaus- 
sian has been used but it is not a good approximation. 
A Lorentzian has also been used but it does not approxi¬ 
mate the actual curve as well as an exponential near the 
origin. Both the exponential and the Lorentzian do well 
at larger lags. We did not use a more complex model 
because the errors caused by weakness in the model are 
dwarfed by the actual variations. 

We fit the auto-covariance to the theoretical model, 
parameterized by the time and frequency scales, tq and 
z/Q and the variances A (white noise) and B (scintillation). 

0(0,0) =A + B (12) 

C(n',0) = Bexp{-\T(n')/Tof^^)(N - n')/N (13) 
0(0, m') = B exjp(—\iy(m')/iyo\)(M — m')/M (14) 

For all the PPTA observations the time scale tq ^ I min, 
so the white noise delta function (A) is clearly separated 
from the temporal scintillation. However the bandwidth 
Vo is often near the frequency resolution (channel width) 
of the receiver and the white noise delta function is not 
always clearly separated from the scintillation. So we 
fit for A, B and tq in the temporal covariance, equation 
(12), first, then we fit for vq in the frequency covariance, 
equation (13), holding A and B fixed. 

4. OBSERVATIONS 

The fluctuations in DM, SDM(t), were determined 
with the technique discussed in Keith et al (2013). The 
fluctuations can be measured with a precision of the or¬ 
der of 1:10®, which is much higher than the absolute ac¬ 
curacy with which DM can be measured. The problem 
with absolute DM is that the pulse shape changes with 
frequency and these changes make it impossible to com¬ 
pare the group delay at different frequencies precisely. 
However the pulse shape is remarkably stable with time. 
Thus we can measure very small time variations even 
though we do not know the absolute DM with compara¬ 
ble accuracy. Here we show SDM (t) with respect to zero 
and provide the mean DM separately. 

The distance of the pulsar is not known precisely in 
either of the two ESE observations. This uncertainty 
propagates into the velocity because, although the an¬ 
gular proper motion is well determined from the timing 
model, it must be multiplied by the distance to obtain 
Vp. In both cases the distance computed from the Tay¬ 
lor and Cordes (1993) model differs significantly from 
that computed using the newer Cordes and Lazio (2001) 
model. Accordingly we show the scattering models de¬ 
rived from both distance estimates to illustrate the sen¬ 
sitivity of the analysis to distance. 
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Figure 2. Observations of J1603—7202. The ESE between 53740 
and 54000 is marked by boxes indicating the parameters used in 
the analysis. Two vertical lines align the scattering observations 
with an obvious step in the DM(t). 

4.1. ESEs 

One of the ESEs, in observations of J1603—7202, is in 
the PPTA data release 1 (Manchester et al. 2013) and 
the dispersion variations were discussed by Keith et al. 
(2013). The other ESE, in observations of J1017—7156, 
was found in the HTRU survey and the dispersion varia¬ 
tions were noted by Ng et al, (2014). It was included in 
continuing PPTA observations, as yet unpublished, but 
taken with the same instruments and the same primary 
analysis procedures as discussed in these earlier papers. 

The observations for J1603-7202 are shown in Figure 
2. The ESE between MJDs 53740 and 54000 is clear 
in all three parameters, 6DM, tq, and z/q, but best de¬ 
fined in vq. The ESE parameters, shown as a solid box 
in the graph, were fit by eye because neither the obser¬ 
vations, nor the analysis are precise enough to justify 
a model fit. The pulsar parameters, the ESE measure¬ 
ments, and the derived parameters are given in Table 
1 labeled (ese). Ignoring Ve we find 0.39 < C < 0.46, 
i.e. the ESE is roughly midway between the pulsar and 
the Earth. If Sout = Wt, we find Wz = 29Wt, and 
3.4 < (rie) < 4.1 cm“^. This lends support to shell 
type models including the corrugated reconnection sheet 
of Pen and Levin (2014). It is also possible that the ESE 
is spherical [Wt = Wz = Sout) but not fully turbulent 
with rms electron density Ri 20% of the mean density. In 
this case (rie) Ri 99 cm“^. This rie is so high that we 
searched for other observations, such as Ha, that might 
show it. We did not find anything interesting, but the 
ESE is so small that it might require a targeted observa¬ 
tion. At this time we regard the spherical ESE model as 
less likely than the extended model. 

We have applied equation 9 to the full line of sight, 
exclusive of the ESE. The results are in Table 1 labeled 
(los). Ignoring Ve we find 0.38 < ^ < 0-45 and the 
resulting average outer scale for the entire line of sight 
is Sout = 1-0 pc. This is consistent with the work of 
Haverkorn et al, (2004, 2008). 


The variation in tq is partially due to significant 
changes in t4ff caused by a combination of Ve and the 
orbital velocity of the binary system. The variation in 
vq is not affected by the velocity and must be due to real 
variations in the turbulence level of the ISM. These are 
very substantial variations and add to the accumulating 
evidence that the ISM is far from homogeneous on AU 
scales. 

The observations for J1017—7156 are shown in Fig¬ 
ure 3. Here the ESE between MJDs 55650 and 55800 
is less well-defined than the one in J1603—7202, but the 
recovery step after the ESE is very abrupt in tq and 
z/Q. The parameters of the analysis are given in Table 1. 
In this case we are inclined to put more weight on the 
Cordes and Lazio (2002) Galactic model because it is 
newer and the Vp calculated with the Taylor and Cordes 
(1993) model is exceptionally high. Ignoring Ve we hnd 
C < 0.17, i.e. the ESE is closer to the pulsar than the 
Earth. If Sout = Wt, we find Wz = 6Wt, so (rie) = 3.7 
cm“^. Again a shell or corrugated reconnection sheet 
would be favored, but a spherical model is tenable if the 
rms electron density is less than the mean by a factor 
Ri y/5. In this case (rie) = 22 cm“^. 

We have applied equation 9 to the full line of sight, 
exclusive of the ESE. The results are shown in Table 1. 
Ignoring Ve we hnd s = 0.39 and the resulting “effective” 
outer scale for the entire line of sight is Sout = 12 pc. This 
is somewhat larger than expected but not inconsistent 
with the work of Haverkorn et al, (2004, 2008). 

One can see that pq is less variable in this pulsar. 
The bandwidth is slightly smaller which provides more 
degrees of freedom in the covariance estimate, but not 
enough to explain the signiHcant reduction in variability. 
This must be a difference in the ISM itself. There is also 
an indication of a very short ESE near MJD 56660. Since 
this possible event is dehned only by a couple of samples 
we will not analyze it further. 

These calculations are only accurate to hrst order. In 
particular the parameters of the ESE are approximate 
and the DM distance is quite uncertain because of un¬ 
certainties in the Galactic rie models. The calculations 
could be improved by including Ve, the binary veloc¬ 
ity of the pulsar, the local velocity of the ISM, and the 
anisotropy of the turbulence as was done for the double 
pulsar J0737—3039A/B (Rickett et al, 2013). We hope 
to do such an analysis for several PPTA pulsars, but it is 
not necessary to obtain a first order model of the ESEs 
reported here. Our goal is to establish that the DM and 
diffractive scattering observations are consistent with a 
simple model of the ESE and provide a good estimate of 
the outer scale of the turbulence in the ESE. 

4.2. Other Events 

One can see an abrupt step in DM(t) for J1603—7202 
in Figure 2 at MJD = 56200, which not clearly present 
in Tq or in eq. We have double-checked the observations 
and see no reason to believe that this step is spurious. 
Both To and eq are extremely variable so perhaps a cor¬ 
relation is simply buried in the natural variation. Never¬ 
theless the integrated electron density changes by -0.002 
pc cm“^ in less than 75 days or 2.6 AU. Glearly the local 
density change must be of the same order as that of the 
earlier ESE. Since the density decreases the earlier idem 
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Table 1 

Model parameters for ESEs. 
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pc 
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days 
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SO 
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Wt 

AU 

Wz 

AU 

tgrp 

Sout 

{rie) 

_‘a 

cm 

J1603-7202(ese) 

1640^ 

61 

38 

0.0023 

260 

0.6 

5 

0.46 

9805 

4.9 

142 

5.1/is 

4.9AU 

3.4 


1170^ 

44 






0.39 

8089 

4.0 

117 


4.0AU 

4.1 

J1603-7202(los) 

1640 

61 

38 

- 

- 

10 

20 

0.45 

39960 

- 

- 

84ms 

l.Opc 

- 


1170 

44 






0.38 

32684 

- 

- 


l.Opc 

- 

J1017-7156(ese) 

30002 

144 

94 

0.0015 

200 

2.0 

2.5 

0.17 

18030 

13.9 

83 

3.3/rs 

13.9AU 

3.7 


8000^ 

384 






0.35 

37698 

29 

174 


29AU 

1.8 

J1017-7156(los) 

3000 

144 

94 

- 

- 

10 

10 

0.39 

52892 

- 

- 

208ms 

12pc 

- 


8000 

384 


- 

- 



0.63 

85671 

- 

- 


9.1pc 
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DM* in pc cm by Taylor and Cordes (1993) DM model; by Cordes and Lazio (2002) DM model. 
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Figure 3. Observations of J1017—7156. The ESE between 55600 
and 55800 is analyzed. A smaller event at 56640 appears to be real 
but it not analyzed in detail. 

tified ESE must be part of a more extended high density 
“local cloud”. Although there are no other statistically 
significant changes in DM(t) there are significant changes 
in the bandwidth, which can only be due to changes in 
the turbulence level. 

A “pinhole” in the ISM can be seen in observations of 
J1713+0747 shown in Figure 4. This pinhole was evident 
in the observations of Keith et al. (2013) but was not 
discussed because it consists of a single DM(t) estimate. 
However it must be real because it is also seen in obser¬ 
vations from Arecibo, Greenbank and Nancay which will 
be discussed in a paper in preparation (Lentati, private 
communication). An under-dense region corresponds to 
a converging “lens”, so one might expect some focussing. 
Accordingly we have also displayed the flux density on 
Figure 4. The problem of seeing small flux changes in 
the presence of strong diffractive scintillation is partic¬ 
ularly clear here. However there is weak evidence for 
a very short jump in i/g a-nd tq at the DM(t) drop of 
0.01 pc cm“^, which would indicate that lower scatter¬ 
ing accompanied the lower density. Of course density 
cannot go negative, so such a small pinhole must be an 
ESE with a hole in it. The observations in fact show a 
marginally significant {2a) jump, followed by the more 
significant drop and recovery to the pre-event level. A 
shell expanding across the line of sight faster than the 
Alfven speed would have such an effect. Such events are 


Figure 4. Observations of J1713-I-0747. The pinhole discussed 
is marked by vertical lines centered on 54800. 

very common in the solar wind, they push up a density 
compression in front of the shock and leave a rarefaction 
behind. 

Nine of the 20 PPTA pulsars showed a clear linear 
gradient in DM(t) in the Keith et al (2013) work. How¬ 
ever in the 3.5 yrs since the end of PPTA DRl, five 
of these pulsars have either lost their gradient or re¬ 
versed it. An interesting example is in the observations 
of J1939-I-2134 shown in Figure 5. Here a steep negative 
gradient changed to an even steeper positive gradient in 
less than 75 days. Remarkably, the negative gradient had 
been observed much earlier (Ramachandran et al, 2006) 
and had remained sensibly constant for 20 years before 
the PPTA DRl. Correlations between all three param¬ 
eters are obvious to the eye. Some of the variation in 
DM(t) is caused by 14; (Keith et al, 2013). As Vp = 13.6 
km/s is unusually low it is likely that Vjsm is also impor¬ 
tant. Although J1939-I-2134 has been analyzed in detail 
(Ramachandran et al, 2006), it would appear that fur¬ 
ther analysis, including the dynamic spectra, as was done 
for the double pulsar by Rickett et al, (2014) would be 
useful. The correlation between pq and tq suggests that 
anisotropy and perhaps changes in Vjsm should be con¬ 
sidered. Indeed the many “bumps” on the DM(t) mea¬ 
surements for J1939-I-2134 may be independent clouds, 
each with different distance, density, turbulence, and ve¬ 
locity. 

Of the PPTA pulsars, the only one which presently 
shows a very linear DM(t) gradient is J1909—3744. The 
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Figure 5. DM(t) observations of J1939+2134. The gradient 
reversal is marked by a vertical line at 55600. 


gradient for this source could be explained by the radial 
velocity of the source as proposed by Cordes (2015). The 
radial velocity of the pulsar’s white dwarf companion has 
been measured at -48±15 km/s (Kerkwijk, 2014). A lo¬ 
cal electron density in the vicinity of the pulsar of 5.6 
cm“^ extending over at least 100 AU would be required 
to explain the observed DM(t) gradient of -7.55 x 10“^pc 
cm“^/day. 


5. DISCUSSION 

ESEs were discovered in flux measurements of compact 
extra-galactic sources (AGNs) because of the strong re¬ 
fractive intensity variations. AGNs are too large to show 
diffractive intensity variations so it was unclear if the 
ESEs were turbulent. The implied density of the ESE is 
of order 10s of cm“^, but this density is quite model de¬ 
pendent as it is very difficult to “invert” refractive fluc¬ 
tuations. Subsequently the phenomenon of “parabolic 
arcs” was discovered in pulsar observations. These are 
diffractive scintillations driven by small scale turbulence. 
These parabolic arc observations show many examples of 
discrete scattering structures of comparable size and rms 
density to ESEs, often many such examples in a single 
observation. It is tempting to identify both types of ob¬ 
servation as ESEs, but those found in parabolic arcs are 
significantly more common and the observations do not 
provide an estimate of the mean density on larger scales. 

Here we show two observations of fluctuations in DM 
accompanied by strong diffractive scintillations, which 
may provide the missing link. The DM(t) observations 
provide a direct measurement of the density, less model- 
dependent than refractive flux variations. The observa¬ 
tions do not have sufficient signal-to-noise ratio to show 
parabolic arcs, but they show strong diffractive intensity 
variations. The observations do not show refractive in¬ 
tensity variations because any such variation would be 
obscured by the strong diffractive intensity variations. 
Modeling the diffractive scintillation shows that these 
structures, and presumably all ESEs, are very efficient at 
diffractive scattering because they have an “outer scale” 
of the order of AU, so they are fully turbulent on a much 
smaller scale than the average ISM (which has an outer 
scale of order pc). 

Ghanges in flux, group delay, or DM(t) would be ob¬ 
served if an ESE passes within Og of the line of sight to 
a compact radio source. However parabolic arcs from an 


ESE are visible when the ESE is within r; 50o from the 
pulsar (Gordes et ah, 2006). Thus a snapshot parabolic 
arc observation is of order 25 times more likely to hnd 
an ESE, than is a comparable observation of refraction, 
DM(t), or group delay. 

Several of the PPTA pulsars show evidence of extended 
small scale variations in DM(t) persisting, in the case of 
J1939-I-2134, for 10 years or 100s of AU. It would be in¬ 
teresting to know more about such high density clouds. 
If molecular lines, or HI absorption for example, could be 
observed one might search for velocity structure. ESEs 
could be dynamic structures formed at the boundary of 
clouds. It would also be interesting to know if they show 
other molecular lines suggestive of partial ionization for 
example, as that would alter the nature of the turbu¬ 
lence. If they are related to corrugated reconnection 
sheets, one could search for other evidence for such re¬ 
connection sheets. In one case (Hill et al, 2005) a group 
of four ESEs moved together across the line of sight to a 
pulsar over the course of 26 days, apparently in a linear 
array. In another case (Brisken, et al, 2010) an image 
of an ESE was found to be roughly linear but well offset 
from the line of sight. These cases would be consistent 
with a sheet or a rope topology. 


The Parkes radio telescope is part of the Australia 
Telescope, which is funded by the Commonwealth Gov¬ 
ernment for operation as a National Facility managed by 
CSIRO. S.O. is supported by the Alexander von Hum¬ 
boldt Foundation. 
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